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Quantum computation provides great speedup over its classical counterpart for certain problems. 
One of the key challenges for quantum computation is to realize precise control of the quantum 
system in the presence of noise. Control of the spin-qubits in solids with the accuracy required by 
fault tolerant quantum computation under ambient conditions remains elusive. Here, we quantita¬ 
tively characterize the source of noise during quantum gate operation and demonstrate strategies to 
suppress the effect of these. A universal set of logic gates in a nitrogen-vacancy center in diamond 
are reported with an average single-qubit gate fidelity of 0.999952 and two-qubit gate fidelity of 
0.992. These high control fidelities have been achieved at room temperature in naturally abundant 
diamond via composite pulses and an optimised control method. 


Quantum computations promise solutions of certain 
intractable problems in classical computations, such as 
quantum simulations prime factoring mm, and 

machine learning mm- Recently, exciting progress to¬ 
wards spin-based quantum computation has been made 
with nitrogen-vacancy (NV) centers in diamond, which 
is a promising candidate for quantum computation un¬ 
der ambient conditions [7]. Long coherence time of NV 
centers has been achieved with dynamical decoupling 
technique [8]. Robust single-qubit and two-qubit gates 
have been accomplished ISHU]. NV-NV entanglement 
has been realized showing the scalability of this quantum 
system [13 HI]. Quantum algorithms [niisi and quan¬ 
tum error correction [13 Hi] have been recently reported 
in NV centers. Further improvement of NV center to¬ 
wards realistic quantum computation would require high 
fidelity quantum gates with errors below fault-tolerant 
threshold, which is proposed to be between 10“^ and 
10“^ depending on the noise model and the compu¬ 
tational overhead for realizing quantum gates [TBfl20j . 
While fault-tolerant control fidelity has been reported 
very recently in superconducting qubits [21] , trapped ions 
[73 - and phosphorus doped in silicons [23, it is still of 
great challenge to achieve fault-tolerant fidelity under 
ambient conditions, which is the case for NV center in 
naturally abundant diamond at room temperature. 
It is because that the noise, not only stemming from in¬ 
teractions between the quantum system and the environ¬ 
ment but also induced by imperfect manipulations, limits 
the fidelities of quantum gates. 

In the following, a universal set of high-fidelity quan¬ 
tum gates at the threshold for fault-tolerant quantum 
computations in the NV center system are realized. A 
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novel composite pulse technique has been developed to 
suppress the noises during the single-qubit gates. We 
adopt the quantum optimal control method for design¬ 
ing the two-qubit gate to suppress the effects induced by 
both the spin bath and the imperfect control field. A nov- 
elly designed coplanar waveguide (CPW) with 15 GHz 
bandwidth has been fabricated to minimize the effect 
due to the finite-bandwidth for the microwave pulses. 
Pulse fixing technique is utilized to correct the effect of 
the imperfect generation of the microwave pulses. The 
average gate fidelity of single qubit is measured to be 
0.999952(6) and the fidelity of the two-qubit controlled- 
NOT (GNOT) gate reaches 0.9920(1). Thus we have suc¬ 
cessfully demonstrated a universal set of quantum gates 
with the fault-tolerant control fidelity in diamonds. 

Results 

NV center in diamond. Fig. la depicts the NV cen¬ 
ter in diamond and the energy level structure. The NV 
center consists of a substitutional nitrogen atom with an 
adjacent vacancy cite (V) in the diamond crystal lattice. 
The ground state of NV center is an electron spin triplet 
state ^A 2 , with three sublevels \ms = 0) and \ms = ±1). 
The intrinsic nitrogen nuclear spin results further split¬ 
ting of the energy levels due to the hyperfine coupling. 
In the experiment presented here, the two-qubit system 
comprises the electron spin qubit and ^'^N nuclear spin 
qubit. Electron (nuclear) spin states \ms = 0) and 
\ms = —1) (|m/) = 0) and \mi = -fl) ) are encoded 
as the electron (nuclear) spin qubit ( Fig. la). When 
a 532 nm green laser pulse is applied to NV center, the 
electron spin will be excited to the state, and then 
fluorescence emission can be measured. The optical tran¬ 
sitions are used to initialize and read out the state of the 
electron spin. The polarization of the electron spin is 
measured to be 0.83(2) (see Supplementary Note 8 and 
Supplementary Eig. 9). A magnetic field of 513 Gauss is 
applied along the NV symmetry axis ([111] crystal axis) 
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Figure 1: Characterization of the noises in experiment, a, Schematic atomic and energy levels of the NV center. Laser 
pulses are used for state initialization and readout. Microwave pulses (MW) are used for coherent control of the electron spin. 
Radiofrequency pulses (RF) are used for manipulating the nuclear spin, b, Result of the Ramsey experiment (inset, pulse 
sequence) for the electron spin. The solid black line in the left panel is fit to the experimental data (red circles) with Po{So) 
, which is the probability density distribution of So depicted in the right panel. The decay time of FID is measured to be 
T 2 = 1.68(3) /rs. c, Results of the nutation experiment (inset, pulse sequence) for the electron spin. The stepped MW pulse 
length is set to be 810 ns. The decay time of the nutation is T 2 = 73(7) /rs. The solid black line in the left panel is fit to the 
experimental results (olive diamonds) with Ti(5i) , the probability density distribution of depicted in the right panel. The 
error bars on the data points are the standard deviations from the mean. 


to split the \ms = ±1) energy levels, and to polarize the 
nitrogen nuclear spin (Supplementary Note 5). 

Realization of high fidelity single-qubit gates. 

The high-fidelity single-qubit gate is demonstrated on 
the electron spin qubit. The state evolution of the elec¬ 
tron spin qubit is governed by the Hamiltonian Hideai = 
27ra;i (cos 4>Sx + sin (pSy) in the rotating frame, where (p is 
the microwave phase and Wi is the Rabi frequency. Un¬ 
fortunately, this is only the ideal case. In practice, the 
Hamiltonian is iJprac = ^ttSoSz + 27r(a;i -I- (5i)[cos(0 -I- 
S(p)Sx -I-sin((/)-h(5(/<)S'y], where Sq, Si, and d<p lead to gate 
error. The i5o mainly comes from the Overhauser field 
(due to the interaction with the nuclear spin bath), the 


magnetic field fluctuation, and the instability of the mi¬ 
crowave frequency. The di partly comes from the control 
field error due to the static fluctuation of the microwave 
power. The distortion of the pulse due to the finite band¬ 
width also contributes to the . The 6(p describes the 
distortion of the microwave phase due to the imperfect 
microwave generation. We take the distortion of the am¬ 
plitude and the phase of the microwave pulse as the sys¬ 
tematic errors, which have unchanged deviations from 
the ideal cases during each experimental scan. The other 
errors can be treated as the fluctuation of the experi¬ 
mental parameters. All of these errors limit the fidelity 
of quantum gates. 
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We first quantitatively characterize the distributions 
of errors due to the fluctuation of the experimental pa¬ 
rameters. Since the timescale of variation of Sq and 
Si is much longer than that of the quantum gates, So 
and (5i can be taken as quasi-static random constants, 
and the practical Hamiltonian is rewritten as i?prac = 
27rSoS2+27r(uJi+Si)Sx- The fluctuations of <5o and (5i lead 
to the decoherence of the qubit during quantum gates 
when averaged over repeated experiments. Fig. lb and 
Ic show the experimental distributions of Jq ^^nd (5i de¬ 
rived by Free Induction Decay (FID) and nutation exper¬ 
iments, respectively. The distribution of So, i.e. Fo(So), 
is obtained from the FID signal via the Ramsey experi¬ 
ment, which is shown in the left panel of Fig. lb. The 
oscillation of the FID signal is due to the detuning of 
the microwave frequency. We assume that So satisfies a 
Gaussian distribution Po(So) = exp(—(5o/2cr^)/(o-\/27r), 
where a stands for the standard deviation of the distribu¬ 
tion. The fitting of the FID data is based on distribution 
Tb(^o)) with parameter a optimized to achieve best agree¬ 
ment between the fitted and experimental data. The fit¬ 
ted data with best agreement is shown as the solid curve 
in Fig. lb, which gives a = 0.131(5) MHz, with the prob¬ 
ability density distribution Po{So) depicted in the right 
panel of Fig. lb. The distribution of (5i, i.e. Pi(^i), 
obtained via nutation experiments, is shown in Fig. Ic. 
The Rabi frequency is set to be uJi = 10 MHz. The Si sat¬ 
isfies a Lorentzian distribution of T’i(^i) = ')/Tr{S'l -l-y^), 
where 7 is the half-width at half-maximum of the dis¬ 
tribution. The distribution Hi((5i) together with PoiSo) 
obtained from the FID, is used to fit the nutation experi¬ 
ment. Best agreement between the fitted and experimen¬ 
tal data is achieved with 7 = 0.0024(4) MHz. The fitting 
result is the solid curve in the left panel of Fig. Ic. The 
probability density distribution Fi((5i) is depicted in the 
right panel of Fig. Ic. 

Since distributions of (5o and (5i have been quantita¬ 
tively characterized, we now demonstrate the realization 
of high-fidelity gates with a novel composite pulse, which 
is designed to suppress these noises simultaneously. To 
form the new composite pulse, we replace one of the 
constituent elementary pulses of CORPSE (compensa¬ 
tion for off-resonance with a pulse sequence) by BBl 
(broadband number 1) sequence. This new composite 
pulse is named as BBlinC for short. Theoretical calcula¬ 
tion shows that BBlinC is robust against both errors 
of do and di (Supplementary Note 3 and Supplemen¬ 
tary Fig. 1). We experimentally compare the perfor¬ 
mance of this new pulse sequence with three other types 
of pulse sequences, which are naive (rectangular) pulse, 
SUPCODE (soft uniaxial positive control for orthogonal 
drift error) [24] pulse, and BBl pulse [25]. The rota¬ 
tion of an angle 9 around the axis in the equatorial plane 
with azimuth (jj is denoted by (0)^. The naive pulse is 
very sensitive to the errors do and di, with leading or¬ 
ders of both errors preserved in the evolution operator 
(corresponding to second orders in gate fidelity). The 
SUPDODE pulse, a type of dynamically corrected gate. 


has been proposed to suppress the dephasing noise dur¬ 
ing quantum gates. The five-piece SUPCODE pulse has 
been used here, where the waiting time D( 3 ) and the 
pulse duration T 2 satisfy the specific requirement (see 
Supplementary Note 3). Under the five-piece SUPCODE 
pulse, up to second order of do can be canceled (corre¬ 
sponding to sixth order preserved in gate fidelity). This 
has recently been demonstrated in NV centers m with 
control field wi of 1 MHz. However, if we increase the 
control field to shorten the gate time, the fluctuation 
of the control field will dominate the decrease of the 
gate fidelity. A composite pulse, named BBl, which is 
shown to be robust against the di, can be applied to 
suppress the fluctuation of the control field. This pulse 
sequence is (9/2)o - (n)^ - ( 271 ) 3 ^ - {n)^ - {9/2)o, with 
(j) = arccos(— 0 / 47 r). The error induced by di is inhib¬ 
ited up to second order in the evolution operator. The 
sequence of BBlinC, which consists of seven pulses, is 
shown in Fig. 2a. The BBlinC sequence is depicted as 
(6»/2)o - (7r)0- (27r)3,^- ( 7 ^) 0 - ( 03 ) 0 -{ 02 )^- {9i)o, where 
(j) = arccos(— 0 / 47 r), 9i = 6*/2 — arcsin[sin( 0 / 2 )/ 2 ], 62 = 
27r—2 arcsin[sin( 0 / 2 )/ 2 ], and 6*3 = 27r—arcsin[sin(0/2)/2]. 
Leading orders of both the do and di errors are canceled 
in the evolution operator, so BBlinC can suppress both 
errors simultaneously. 

The average gate fidelity |2S| between a quantum op¬ 
eration ^ and a target unitary quantum gate U is defined 
as 


F.itu)= / 


( 1 ) 


where the integral is over the uniform measure dip on the 
Hilbert space of the system. To estimate the fidelity by 
equation (14), it is necessary to characterize the quan¬ 


tum operation ^ with comprehensive knowledge about 
the Hamiltonian of the quantum system, and the Hamil¬ 
tonian of the control together with the errors during the 
gate operation. We numerically calculate the gate fideli¬ 
ties (Supplementary Note 2) for these pulse sequences 
with both do and di . The BBlinC pulse presents a de¬ 
sirable region with fidelity higher than 0.9999 (Supple¬ 
mentary Note 3 and Supplementary Fig. 1). The control 
field strength is set to uji = 10 MHz for experimental 
realization of high-fidelity single-qubit gates, so that the 
error due to do and di can be greatly reduced. 

Though composite pulses provide possibility to realize 
high fidelity gates in solid-state system, the performance 
of the composite pulses is limited by the systematic er¬ 
rors, which are usually not taken into account in the 
composite pulse designing. Herein, we presents methods 
to correct the systematic errors. The distortion of the 
microwave phase is corrected by pulse fixing technique, 
which is previously developed in liquid NMR (nuclear 
magnetic resonance ) m- The phase of the microwave 
pulses was recorded by an oscilloscope, and then the de¬ 
viation of the phase from the ideal case was fed back 
to the arbitrary waveform generator so that this distor¬ 
tion can be minimized (see Supplementary Note 6 and 
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Figure 2: Characterization of the performance of single-qubit gates, a, Pulse sequences corresponding to naive 
pulse, SUPCODE pulse, BBl pulse and BBlinC pulse. The time durations (i = 1, 2 and 3) of SUPCODE is included in 
Supplementary Note 3. b, Results of randomized benchmarking. Red diamonds, blue circles, navy blue squares, and green 
triangles represent measured single-qubit hdelity of the output state from each of the individual sequences of gates. Solid lines 
are hts of equation ([^ with averaging each over all the randomization pulse sequences. The average gate fidelity for the naive, 
five-piece SUPCODE, BBl and BBlinC pulses are 0.99968(6), 0.99916(8), 0.999945(6), and 0.999952(6) respectively. The inset 
shows the average error per gates (sg) of the pulses. The error bars on the data points are the standard deviations from the 
mean, and those of £g in inset are given by errors of the ht of the randomized benchmarking data. 


Supplementary Fig. 5). The distortion of the ampli¬ 
tude is primarily due to the limited bandwidth of the 
microwave fed structure. We designed and fabricated an 
ultra-broadband CPW with a bandwidth up to 15 GHz 
(see Supplementary Note 5 and Supplementary Fig. 4). 
Then the effect of the finite bandwidth is found to be 
negligible. The reflection between the microwave com¬ 
ponents also contributes to the imperfection of the mi¬ 
crowave pulses. This effect was diminished by inserting 
proper attenuators between the microwave components 
(Supplementary Note 6 and Supplementary Fig. 6). The 
amplitude distortion is further corrected with pulse fix¬ 
ing (Supplementary Note 6 and Supplementary Fig. 7). 
After all these procedures, the microwave pulses can be 
fed to the electron spin with almost perfect pulse shapes 
. The detailed procedures for correcting the systematic 


errors are included in the Supplementary Note 6. 

To experimentally quantify the performance of the 
single-qubit gates, we utilize the randomized bechmark- 
ing (RB) method [28]. In the RB experiment, quantum 
gates are evaluated by measuring the performance when 
random sequences of the gates are applied. The error 
due to the imperfect measurement and state preparation 
can be separated, and the gate fidelity can be determined 
precisely. When the number of randomized gates imple¬ 
mented in a sequence is increased, the accumulated gate 
error reduces the measured single-qubit fidelity F of the 
output state, which is defined as the overlap of the ideal 
and the measured states. The average output state fi¬ 
delity F |2S] can be written as 

F = 1/2 +1/2(1-dif)(l-2£g)( 


( 2 ) 
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Figure 3: Optimization of the control of the two-qubit system to achieve high fidelity control in the presence 
of noises, a and b, Comparison of the original GRAPE (a) with the modified one (b) in the presence of the noises 5o and . 
Here stands for the maximal strength of the control field in the pulse sequence, and n stands for the number of iterations 
in the algorithms. The upper figures of (a) and (b) are the fidelity contour maps of the two GRAPE methods. The area with 
a control fidelity higher than 0.99, is larger for the modified method than the original one. c. The schematic diagram of the 
optimal control pulses used in the experiment, with the left and right being the amplitude and phase sequences, respectively. 


where da is due to the imperfection of the state initializa¬ 
tion and readout, Eg is the average error per gate, and I is 
the number of randomized quantum gates. The average 
gate fidelity (Tkvg = 1 — Eg) derived from this method 
is resilient to the state preparation and measurement er¬ 
rors. The detailed procedure of RB is included in the 
Supplementary Note 9. Fig. 2b shows the result of the 
randomized bechmarking, where the control field wi is set 
to be 10 MHz. The average gate fidelity of naive pulses is 
found to be 0.99968(6). The five-piece SUPCODE pulses 
provide an average gate fidelity of 0.99916(8) which is 
lower than that of the naive one. This is because the fluc¬ 
tuation of control field dominates the decay of the fidelity 
when wi = 10 MHz. The average gate fidelity is greatly 
improved when BBl pulses are applied. The extracted 
average gate fidelity for BBl pulses is 0.999945(6). The 
highest fidelity is achieved when BBlinC pulses are used. 
The average gate fidelity for BBlinC is 0.999952(6), pro¬ 
viding an error per gate below lO”'^ which is about one 
order of magnitude lower than that of naive pulses (Sup¬ 
plementary Table 1). 

Realization of high fidelity CNOT gate. Two- 


qubit CNOT gate together with the single-qubit gates 
provide a universal set of quantum gates. Two-qubit 
gates have been demonstrated in NV centers mm- 
However, realization of the two-qubit gates with high fi¬ 
delity to meet the requirement of the fault-tolerant quan¬ 
tum computation is still of great challenge. The control 
qubit is the nuclear spin qubit and the target qubit is 
the electron spin qubit. The detailed Hamiltonian is in¬ 
cluded in Supplementary Note 1. The CNOT gate is 
designed to flip the electron spin qubit iff the nuclear 
spin qubit is in the state |m/ = 1). Microwave (MW) 
and radio-frequency (RF) pulses are used to manipulate 
this two-qubit system. 

To achieve high fidelity, we have improved a type of 
quantum optimal control method, named GRAPE (gra¬ 
dient ascent pulse engineering) [29], to design the CNOT 
gate which is robust against both Sq and 5i. GRAPE 
has recently been used to realize the quantum error cor¬ 
rection |16j and a high-fidelity entanglement m in dia¬ 
mond. Because the nuclear spin is insensitive to the noise 
(such as fluctuations of the external magnetic field and 
the control RF field) and the CNOT gate designed by the 
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quantum optimal control method consists of microwave 
pulses only, we can take into account the noise felt by the 
electron spin in the optimization procedures. We modi¬ 
fied the GRAPE algorithm to design target gates which 
are robust against both the errors So and Si (see Sup¬ 
plementary Note 4 and Supplementary Fig. 2). Fig. 3 
shows the optimization of microwave pulse parameters 
to achieve the high control fidelity. Fig. 3a (b) shows 
the optimizing procedure of original (modified) GRAPE 
algorithm. The final achieved high fidelity regions with 
the two GRAPE algorithms are compared in the upper 
parts of Fig. 3a and 3b. It is clear that the modified 
method presents much more robustness against both i5i 
and So than the original one. The left and right pan¬ 
els of Fig. 3c show the amplitude and phase sequences 
after the improved optimization, respectively. The to¬ 
tal length of the sequence for one GNOT gate is 696 ns, 
which consists of twelve pulses of equal length with dif¬ 
ferent amplitudes and phases. The theoretical fidelity of 
the GNOT gate via this pulse sequence is estimated to 
be 0.9995 in the absence of So and di. If the noises i5o 
and i5i are taken into account, the average gate fidelity 
of the designed GNOT gate is estimated to be 0.9927. 
The detailed pulse optimization and fidelity estimation 
are included in the Supplementary Note 4. 

Fig. 4 presents the experimental realization of the 
GNOT gate. The two-qubit system is prepared to the 
state |01) and |00) in Fig. 4a and 4c, respectively (The 
state \ms = 0,m/ = 1) is denoted as |01) and the state 
\ms = 0,TO/ = 0) as |00) hereafter). The black (red) 
symbols with error bars in Fig. 4a and 4c are the exper¬ 
imental FID signal via the Ramsey experiment without 
(with) the GNOT gate. Fig. 4b and 4d are the fast 
Fourier transformation of the FID signals, so that the re¬ 
sult of the GNOT gate can be observed in the frequency 
domain. The two peaks correspond to the nuclear spin 
qubit states \mi = 0) and |m/ = 1). The distance be¬ 
tween the two peaks is due to the hyperfine coupling. 
It is clear that when the nuclear spin qubit is of state 
\mi = 1), the electron spin qubit is rotated by tt ( Fig. 
4a and 4b). If the nuclear spin qubit is of state |m/ = 0), 
the electron spin qubit remains unchanged (see Fig. 4c 
and 4d). 

Two-qubit randomized benchmarking, which can be 
used to measure the fidelity of two-qubit gate, requires 
operations on both qubits m- In the hybrid system com¬ 
posed of electron and nuclear spins, single-qubit gates 
of the nuclear spin cost much longer time than that of 
the electron spin qubit. The typical operation time on 
the nuclear spin qubit (about 50 ^s for a tt rotation) 
is much longer than the dephasing time of electron 
spin qubit. The decoherence effect on the electron spin 
during the nuclear spin operation in the two-qubit RB 
experiments will dominate the fidelity decay in random¬ 
ized benchmarking, and the gate fidelity of GNOT can 
not be precisely determined in this way. On the other 
hand, though process tomography provides a full charac¬ 
terization of GNOT gate, the measured gate fidelity with 


this method is sensitive to errors in state preparation 
and measurement. In this scenario, repeated application 
of GNOT gates on the system and recording the dynam¬ 
ics of the quantum state [SD] will be a good choice to 
estimate the gate fidelity. 

Fig. 4e shows the population of |01), P|oi) , after re¬ 
peated applications of GNOT gates on the input state, 
which is generated by a selective tt/2 RF pulse. The inset 
of Fig. 4e shows the pulse sequence, where the number 
of repeated GNOT, N, is even. When N is increased, 
the error of the GNOT gate will accumulate and T’loi) 
will decay. A wealth of information can be obtained by 
studying the state dynamics under repeated applications 
of GNOT gates. In ref. |30l the dynamics of states obey 
a simple exponential decay and the gate fidelity can be 
extracted with twelve applications of GNOT gates. How¬ 
ever, when more GNOT gates (up to 192 in this work) 
are applied, we find that the decay is not exponential. 
So the gate fidelity cannot be simply extracted from the 
decay as presented in the ref. [301 In Fig. 4e, the dynam¬ 
ics of the state exhibits oscillatory attenuation as N in¬ 
creases. The oscillation results from the deviation of the 
realistic evolution from the ideal GNOT operation. For 
example, the optimization procedure adopts the hyper¬ 
fine coupling A = —2.16 MHz, which may differ slightly 
from the experimental value Aexp by SA = A^^p — A. The 
frequency of the microwave will not be equal to the res¬ 
onant frequency exactly, which induces an off-resonance 
term (JH) in the practical Hamiltonian. The decay of 
the T’loi) is mainly due to the static fluctuation of di and 
(5o, which can be quantitatively characterized by the FID 
and the nutation experiment, respectively. The values of 
SA and S^l were derived to be 0.008 MHz and 0.068 MHz, 
by fitting the experimental data (Supplementary Note 9 
and Supplementary Fig. 10). The fitting result shown 
as the solid blue line in Fig. 4e agrees with the exper¬ 
imental data. Since the comprehensive information on 
the Hamiltonian of the quantum system and the control 
field together with the errors is ready, the gate fidelity 
of GNOT can be directly obtained according to equation 
(14) with 0.9920(1) (see Supplementary Note 9). 

Discussion 


To sum up, we have achieved fault-tolerant control fi¬ 
delity for a universal set of quantum gates in diamond 
with of natural abundance under ambient conditions. 
Several errors which limit the fidelity of quantum gates 
have been quantitatively characterized and effectively 
suppressed. With new composite pulses, we have realized 
single-qubit gate with fidelity up to 0.999952. A modi¬ 
fied optimal control method has been developed to design 
the control pulse for the GNOT gate, which achieves a 
gate fidelity of 0.992. To the best of our knowledge, our 
results stands for the state of art high-fidelity control of 
solid-state spins under ambient conditions. Our method 
can not only be used to realize high fidelity GNOT gate 
in the system consisting of on the electronic spin and 
the host nitrogen nuclear spin in NV, but also can be 
applied for the coupled NV-NV systems, which provide 
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Figure 4: Characterization of the performance of CNOT gate, a and c are experimental data and the fittings of FID 
(black : without applying the CNOT gate, red: with applying the CNOT gate), when the initial states are |01) and |00), 
respectively. Lines are fits to the experimental data (symbols with error bars), b and d are the fast Fourier transform results 
of FID signals, e. The probability in the state |01) (T'|oi>) after applying CNOT gates repeatedly, with the sequence shown 
in the inset. The fitting (blue line) agrees with the experimental data (blue circles with error bars). The error bars on the 
data points are the standard deviations from the mean. The detailed procedure of normalization to obtain P|oi> is included in 
Supplementary Note 7 and Supplementary Fig. 8. 


possibilities for future scalable architectures. The noises 
which limit the fidelity of the two-qubit gate in NV-NV 
system can also be greatly suppressed and high-hdelity 
gates are available by our method (Supplementary Note 
10 and Supplementary Fig. 11). The methods presented 
here to achieve high control fidelity can also be applied 
to other quantum systems, such as quantum dots, phos¬ 
phorus doped in silicon, and trapped ions. 

Methods 

Experiment setup. The NV center in [100] face bulk 
diamond was mounted on a typical optically detected 
magnetic resonance (ODMR) confocal setup which was 
synchronized with the microwave bridge by a multichan¬ 
nel pulse generator (Spincore, PBESR-PRO-500). The 
nitrogen concentration was less than 5 ppb and the abun¬ 
dance of was at the natural level of 1.1%. The 532 
nm green laser for pumping and phonon sideband fluo¬ 


rescence ( 650-800 nm) went through the same oil objec¬ 
tive (Olympus, PLAPON 60X0, NA 1.42). In order to 
preserve the NV center’s longitudinal relaxation time Ti 
from laser leakage effects, the pump beam was passed 
twice through an acousto-optic modulator (ISOMET, 
power leakage ratio ~1/1000) before it went into the ob¬ 
jective. We created a solid immersion lens (SIL) in the 
diamond around an NV center (see Supplementary Note 
5 and Supplementary Fig. 3). The SIL increases the PL 
rate to about 400 kcounts s“^. The fluorescence inten¬ 
sity was collected by an avalanche photodiodes (Perkin 
Elmer, SPCM-AQRH-14) with a counter card. The ad¬ 
justable external magnetic field, created by the perma¬ 
nent magnets, was aligned by monitoring the variation 
of fluorescence counts. Microwave and radio frequency 
pulses were carried by ultra-broadband coplanar waveg¬ 
uide with 15 GHz bandwidth. 
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Supplementary Figure 5: Pulse sequences and contour plot of average gate fidelity of single-qubit gate. The upper 
panels show the pulse sequences of ( 6^)0 with (a) rectangular naive, (b) five-piece SUPCODE, (c) BBl, and (d) BBlinC pulses, 
where the phase and duration of each piece of pulse is depicted. The lower panels show the average gate fidelity of (7r/2)o with 
the errors tio and when (a) rectangular naive, (b) five-piece SUPCODE, (c) BBl, and (d) BBlinC pulses are applied. The 
regions of fidelity larger than 0.9999 are surrounded with black curves for clarity. 



Supplementary Figure 6: Designed pulse sequence for CNOT and average gate fidelity of the seqnence. (a) 
Amplitude (left panel) and phase (right panel) of the designed pulse sequence, (b) Calculated average gate fidelity of the 
sequence under the errors (5o and (5i. 
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Supplementary Figure 7: Image of the SIL in diamond. 
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Supplementary Figure 8: Scattering parameters of the CPW. 
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Supplementary Figure 9: Correction of microwave pulse distortions, (a) Schematic correction diagram showing the 
main instruments generating, processing, and sampling the microwave pulses, (b) Extracted phase distortions without (left 
penal) and with (right penal) the correction, (c) Calculated average gate fidelity of a BBlinC 7r/2 gate with sampled waveforms 
without (left penal) and with (right penal) the correction. The region of fidelity larger than 0.9999 is labeled. 
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time (ns) time (ns) 

Supplementary Figure 10: Correction of pulse distortions caused by leakage and reflection, (a) Waveform of a GRAPE 
pulse without improvement. The distortions are mainly caused by leakage and reflection of the diplexer. (b) Improved waveform 
of the pulse sequence. The distortions are corrected by inserting a 10 dB attenuator between microwave components. 
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Supplementary Figure 11: Correction of the amplitude distortions, (a) Waveform without correction, (b) Waveform with 
correction. 
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Supplementary Figure 12: Schematic normalization sequences in the two-qubit experiment. Here pi denotes the 
initialized state with laser, pf denotes the final state after applying control sequence to pi. Trap (ttrpi) is a radio-frequency tt 
pulse driving the nuclear spin transition between states \ms = 0,mj = 1) and \ms = 0,m/ = 0) (transition between states 
\ms = 0,mi = 0) and \ms = 0, mj = ~1))- The measured photoluminescence intensity after each sequence is denoted by Si, 
S 2 , Ri, or R 2 . 
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Supplementary Figure 13: Experimental measurement of the polarization of the NV electron spin. The polarization 
a is extracted from the ratio of the two nutation amplitudes to be 0.83(2). 
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Supplementary Figure 14: Extraction of the practical Hamiltonian. The experimentally measured Ejoi) is shown as the 
black circles. The practical Hamiltonian is extracted by fitting the experimental data. The fitted Eioi) is shown as the red 
solid line. Best-fit values of 5A and 5Q. are 5A = 0.008(1) MHz and tffl = 0.068(1) MHz. For comparison, the blue dashed line 
shows the simulated result with 5A = 0 and = 0. 
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Supplementary Figure 15: Pulse sequence for the CNOT gate in NV-NV coupled system, (a) Amplitudes of the 
control microwave fields, (b) Phases of the control microwave fields. 
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III. SUPPLEMENTARY TABLES 


When MW pulses are applied, the total Hamiltonian be¬ 
comes 


Supplementary Table I: Summarization of the results of 
single-qubit randomized benchmarking. Average gate 
fidelity (Tk) and error per gate (Sg) is shown for naive, five- 
piece SUPCODE, BBl, and BBlinC pulses. 


Pulse sequence 

Fa 

Eg 

naive 

0.99968(6) 

3.2(6) xlO"'^ 

five-piece SUPCODE 

0.99916(8) 

8.4(8) xlO""' 

BBl 

0.999945(6) 

5.5(6)xl0"® 

BBlinC 

0.999952(6) 

4.8(6) X10"® 


H = Hnv + Hc^ 


( 6 ) 


with 


He = 27rv^wi cos(27r/Mwi + (fySx, (7) 


where (f) is the phase of the MW pulse, uji is the amplitude 
of the MW pulse. 

The Hamiltonian can be transformed into the rotating 
frame as 


Hrot — Uti 


sHUl 


iUu 


dULns 

dt 


( 8 ) 


with 


IV. SUPPLEMENTARY NOTES 
A. Supplementary Note 1 

Hamiltonian of the NV system. The NV center 
includes a substitutional nitrogen atom and a vacancy 
in the nearest-neighbor lattice position. In our experi¬ 
ment, a static magnetic field, Bq = 513 G, is applied 
along the NV symmetry axis ([1 1 1] crystal axis). The 
Hamiltonian of the NV center can be written as 


[/trans = (9) 

With rotating-wave approximation, the Hamiltonian in 
the rotating frame can be simplified as 

+ si) - uiih - II) + A{Sl + SJ,)) 

+ He, rot, 

( 10 ) 

where 

He, rot = 2'k[5VLSI + uJi/ V2{cos<j)Hx + sm(j)Hy)], (11) 
with 


Hnv = MDSl + oJsS, + Pll - ioih) + iLhf, (3) 

where uis = —7ei3o/27r (w/ = 7Ni3o/27r) is the Zeeman 
splitting of the electron (^^N nuclear) spin, 7 e ( 7 n) is 
the electronic (^^N nuclear) gyromagnetic ratio, and 
Iz are the electron and nitrogen nuclear spin operators, 
respectively. The zero field splitting D = 2870 MHz and 
the nuclear quadrupolar splitting P = —4.95 MHz. The 
hyperfine interaction between the NV electron spin and 
the ^^N nuclear spin is 

Hi^f = 2Tr[A±{SxIx + Syly) + ASzIz\j (4) 


^0 1 o\ 
10 1 
Vo 1 o) 


H, 


1 

71 


0 -* 0 
i 0 i 
0 -i 0 




( 12 ) 


(13) 


, (5H = D — LOS — A — /mw, /rf = —P + w/, and I 
representing 3x3 identity matrix. 


The strength of the hyperfine interaction is about 2 MHz. 
Because of the strong zero field splitting and Zeeman 
splitting terms of the electron spin, the effect of the in¬ 
teraction term S^Ix + ^yly is strongly suppressed and 
can be neglected. A = —2.16 MHz is determined via 
CW ESR experiment. In the secular approximation, the 
Hamiltonian is 

Hnv = 2tt{DSI + LosSz + ASzh + Pll - w/G), (5) 

The electron (nuclear) spin states \ms = 0) and \ms = 
— 1 ) (|to/ = 0 ) and \mi = - 1 - 1 ) ) are encoded as the 
electron (nuclear) spin qubit. 

Microwave (MW) and radio-frequency (RF) pulses are 
used to manipulate the two-qubit system. The frequency 
of MW and RF pulses are /mw and /rf , respectively. 


B. Supplementary Note 2 

Calculation of average gate fidelity. The average 
gate fidelity between a quantum operation / and a target 
unitary quantum gate U is defined as 

^a(?,G) = J (14) 

where the integral is over the uniform measure dtp on 
state space, normalized so /d')/ = 1 [ 1 ]. 

In the single-qubit case, the average gate fidelity can 
be derived to be [ 2 ] 

= ^ + ^ Y, tr{UayUH{ay)), (15) 

j^x,y,z 
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where ax, ay, and a^ are Pauli matrices. 

Quantum optimal control method [3] is used to design 
the pulse sequence of CNOT gate. To calculate the aver¬ 
age gate fidelity of this CNOT gate, Eqn. [^is general¬ 
ized so that the integral is on the two-qubit space [3] . The 
nuclear spin is much less sensitive to the external mag¬ 
netic noise than the electron spin and the GRAPE pulse 
sequence contains only microwave pulses, so the decoher¬ 
ence during the operation mainly comes from the static 
distributions of So and (5i for the electron spin qubit. 
Then the operation can be expressed as 

iiWW) = 

J dSo J dSiPo{So)Pi{Si)U,eo{So,SiM){i;\u!,^{So,6i), 

(16) 

where C/seq(<5o: ^i) is the calculated two-qubit evolution 
according to the pulse sequence, with the errors do and 
di considered in Hamiltonian. Substituting Eqn. |16| into 
Eqn. [^yields the average gate fidelity between the op¬ 
eration ^ and the target CNOT gate Ccnot, 

Fai ^, Ccnot) = 

J dSo J dS,Po{So)Pii6i){tr{MM^) + MM)\^), 

(17) 

with d = 4 and 

M = C/^NOTt^seq(do,di). (18) 

It can be easily obtained from Eqn. [^that the fidelity 
of operation ^ without the effect of the noise (do and di) 
can be written as 

i"seq = ^(^^[MMMt) + |tr(M)n, (19) 

where the values of do and di are zero. 

C. Supplementary Note 3 

High fidelity single-qubit quantum gates. Con¬ 
sidering a single-qubit gate corresponding to a rotation of 
angle 9 around the x axis on the Bloch sphere, such a gate 
can be realized by evolution under the effective Hamil¬ 
tonian idideai = 27rwin • S, where S = {Sx, Sy, Sz) is the 
spin vector operator of the qubit, n is a three-dimensional 
vector, and the strength uii is a real parameter. The av¬ 
erage gate fidelity is limited by interaction of the qubit 
with environment and fluctuation of the control field. We 
consider the model where the Hamiltonian for rotation 
about the x axis under practical conditions is described 
as iJprac = ‘2'TtSqSz + 2tt{uji + di)[cosdf/f^a; -I- sind(/)S'y]. 
The error do in the Hamiltonian is due to the interaction 
of the qubit with environment, the error di is due to fluc¬ 
tuation of the control field strength and phase error Scj) 
is caused by the imperfect microwave pulse generation. 
Phase error can be efficiently eliminated by pulse fixing 


technique (detailed in Section IV FI and we take it as of 
zero value in this section. We consider the case where 
both do and di vary in a timescale much longer than that 
of quantum gates. In this case do and di are taken as 
quasi-static random constants. 


Supplementary Figurej^ shows the performance of the 
gate by simply applying a naive rectangular pulse. Here 
the gate (7r/2)o (we denote the rotation of an angle 9 
around the axis in the equatorial plane with azimuth (j) 
as ( 0 ) 0 ) is taken as an example. The average gate fidelity 
of (7r/2)o is calculated with respect to different values of 
do and di. The naive pulse is vary sensitive to the errors 
do and di, with leading orders of both errors preserved in 
the evolution operator (corresponding to second orders 
in the average gate fidelity). This corresponds with the 
small region of high average gate fidelity shown in the 
lower panel of Supplementary Figure [^. 

Supplementary Figure shows a type of dynamically 
corrected gate, five-piece SUPCODE [5]. The pulse se¬ 
quence is depicted as ti — {9/2)o — T 3 — (9/2)o — Ti for 
9 e (27r,37r). Here ri = cscd(I — 2cos| -I- cos 9 + 

— 8 cos 1+4cos 6* -I- 0 sin 9) and T3 = —2(riC0s| -I- 

sin I) are durations when control field is off. Under the 
five-piece SUPCODE pulse, up to second order of dp can 
be canceled (corresponding to sixth order preserved in 
the average gate fidelity). In the lower panel of Supple¬ 
mentary Figure [^, the average gate fidelity of (2.57r)o 
(equivalent to ( 7172)0 in the single-qubit case) is shown 
as an example. The region of high average gate fidelity 
is largely extended in the axis of dp, compared with that 
by the naive pulae. 

In Supplementary Figurea type of composite pulse, 
BBl [6], is shown. The pulse sequence is {9/2)o — — 

(27r)30 — (tt)^ — { 9 / 2 ) 0 , with (j) = arccos(—0/47r). Un¬ 
der the BBl pulse, up to second order of di is canceled 
in the evolution operator (corresponding to sixth order 
preserved in the average gate fidelity). The average gate 
fidelity of (7r/2)o with the BBl pulse is shown in the lower 
panel of Supplementary Figurej^. It exhibits a larger re¬ 
gion of high fidelity in the axis of di. Thus suppressing 
the di error by applying the BBl pulse, in combination 
with a proper selection of control field strength wi to sup¬ 
press the do error, can contribute to realization of a high 
fidelity (e.g. 0.9999). 

Supplementary Figure shows a pulse sequence sup¬ 
pressing both the do and di errors simultaneously. The 
sequence is designed by incorporating the BBl pulse 
within the CORPSE pulse, and is named BBlinC here 
for short. There are similar pulse sequences to suppress 
the do and di error simultaneously [ 7 ]. The BBlinC se¬ 
quence is depicted as {9/2)o — (Tr)^ — ( 271)30 — ( 7 r )0 — 
(^* 3)0 - (^* 2)77 - (^* 1 ) 0 , where (j) = arccos(-6'/47r), 6»i = 
9/2 — arcsin(sin(6>/2)/2), 02 = 27r — 2arcsin(sin(0/2)/2), 
and 03 = 27r —arcsin(sin(0/2)/2). Leading orders of both 
the do and di errors are canceled in the evolution opera¬ 
tor. The lower panel of Supplementary Figure shows 
the average gate fidelity of (7r/2)o with the BBlinC pulse. 
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The region of high fidelity is much larger, extended in 
both axes of Sg and i5i. 

Recently the robustness of composite pulse sequences 
against time-dependent noise is analyzed [5]. It is shown 
that composite pulses may also be successfully employed 
in the presence of time-dependent noise. The robustness 
against static as well as time-dependent noise enables 
composite pulse an effective method to improve single¬ 
qubit gate fidelity. 


D. Supplementary Note 4 


Quantum optimal control method for designing 
the CNOT gate. GRAPE is a type of quantum opti¬ 
mal control method. It can be utilized to design control 
sequence to realize a target gate with high fidelity. The 
control sequence contains JV piece of pulses, with the am¬ 
plitude and phase of each piece being different. The total 
Hamiltonian of the fcth pulse in the rotating frame is (see 
Eqn. 10 and Eqn. 11) 


Hrot.fc = 27r(uJs(S, + - J2) + A(S^ + SJ,)) 

T ^C, rot.fci 

( 20 ) 

with 


He , rot,A; = 2TTUJi^k/V2{cos(j)kH^ + sm(j)kHy), (21) 

where (j)k is the phase of the kth pulse, uji^k is the ampli¬ 
tude of the fcth microwave pulse and = 0. 

The evolution operator under the fcth pulse is written 
as 


C/fe = (22) 

where tk is the duration of the fcth pulse. The total 
evolution under the entire sequence is 

1 

Utot =l[Uk. (23) 

k=N 

The time duration of each pulse, tk, is set to be equal 
value r. The two-qubit evolution operator can be de¬ 
scribed as 

1 

Useg{Wi,k,M) = n (24) 

k=N 

where V is the projection operator on the two-qubit sub¬ 
space. 

The target gate is the CNOT gate. 


The performance function of the GRAPE algorithm is the 
fidelity ELeq, which is a function of The values 

of {uJi k, f/'fc} are initialized with random numbers within 
the experimental conditions. The performance function 
are maximized by iteration. Within each iteration, the 
value of Fseq{{uJi,k, 'Pk}) as well as its derivative to uji^k 
and (pk is calculated, then the value of uii^k iPk) is re¬ 
placed by the result of its previous value plus the deriva¬ 
tive multiplied by a proper coefheient. 

This method works well for designing a sequence with 
high fidelity, if errors due to qubit-environment interac¬ 
tion and control field fluctuation are not taken into ac¬ 
count. However, our aim is to realize quantum gates, 
which are not only of high fidelity, but also being robust 
to the errors. A method has been presented [3] to design 
pulse sequence which is robust against the inhomogeneity 
of wi. Herein, we generalize this method to design pulse 
sequence, which is not only robust against the inhomo¬ 
geneity (5i, but also is insensitive to the dephasing noise 
(5o. The performance function of the modified GRAPE is 
defined as Ea(£:, Hcnot) (see Eqn. 17). The new perfor¬ 
mance function is maximized by iteration of the GRAPE 
algorithm. In practical implementation of the modified 
algorithm, the integral in Eqn. [^is replaced by sum of 
discrete points. We find that three points of Jq (^i) are 
enough. 

Supplementary Figure shows the amplitude and 
phase of the designed pulse sequence. The sequence con¬ 
sists of twelve pieces of pulse. The duration of each pulse 
is 58 ns. Without considering the errors Sg and i5i, the 
sequence produces a two-qubit operation Ucai, 


( -0.0060 


^cal — ^ 


1.5974i 


V 


0 1 0 \ 

0 0.9996 0 0.0154 

0.9999 0 0.0059 0 

0 -0.0146 0 0.9991 j 


+ ie 


1.5974j; 


( 0.0032 
0 

0.0119 

V 0 


0 

0 

0 

-0.0253 


-0.0013 

0 

0.0032 

0 


0 \ 
-0.0249 
0 

0.0316 


/ 

(26) 

The fidelity of Hcai is 0.9995. Supplementary Figure 
shows the robustness of the sequence against the errors Jq 
and (5i. When the experimental distributions Pg{Sg) and 
Pi{Si) (which are determined from the experiments) are 
considered, the sequence provides an average gate fidelity 
of 0.9927. 


E. Supplementary Note 5 


Gcnot 


/ 0 0 1 0 \ 
0 10 0 
10 0 0 
V 0 0 0 1 / 


Alignment of the magnetic field. We used the flu¬ 
orescence dependence on the misalignment angle to align 
the magnetic field. According to the literature [S], the 
fluorescence of NV center is sensitive to misalignment 
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angle of the NV axis from a magnetic field Bq when 
the magnitude is approximately 513 G. The difference 
in fluorescence counts is still noticeable even when the 
misalignment angle is only 0.5°. In our experiment, the 
fluorescence count was the same (within counting errors) 
for Bq Ki 0 G and Bq « 513 G. So we estimate the mis¬ 
alignment angle to be within 0.5°. 

Creation of a solid immersion lens. All measure¬ 
ments in our experiment are based on detection of the 
NV photoluminescence. Much of the photoluminescence 
is lost at the diamond surface due to internal reflection. 
The problem can be overcome by creating a solid immer¬ 
sion lens (SIL) [TU]. We created a SIL in the diamond 
around an NV center (Supplementary Figurej^. The SIL 
increases the PL rate to about 400 kcounts s~^. 

Ultra-broadband coplanar waveguide. In the 
experiment, manipulation of qubits is achieved by mi¬ 
crowave (MW) and radio-frequency (RF) pulses, which 
are applied through a coplanar waveguide (CPW). The 
ultra-broadband CPW is designed and fabricated. Sup¬ 
plementary Figure shows scattering parameters of the 
CPW. Up to 15 GHz, the S21 parameter is larger than -3 
dB, and the Sll parameter is about or less than -10 dB. 
Such wide bandwith ensures that there is almost no extra 
distortion of microwave / RF pulses with this CPW. 


F. Supplementary Note 6 

The imperfect devices generate microwave pulses with 
non-ideal amplitudes and phases. The imperfection of 
microwave pulses sent to the NV centers are carefully cor¬ 
rected with pulse fixing technique m- Supplementary 
Figure]^ shows the main instruments for the microwave 
pulse generation. The pulses are generated by an arbi¬ 
trary waveform generator (M8190A, Keysight ), and am¬ 
plified with a power amplifier (ZHL-30W-252-S-I-, Mini 
Circuits). The imperfections in the instruments cause 
distortions of the microwave pulses, which may dramati¬ 
cally decrease hdelity of quantum gates. An oscilloscope 
(DSO-X 92004Q) is used to sample the microwave pulses. 
The pulse distortions are then corrected by predistorting 
the pulses in the right way that the predistortions cancel 
with the distortions by the imperfections of instruments. 

Microwave phase correction. Supplementary Fig¬ 
ure & shows the distortion of the microwave pulse phase 
with / without the correction. It is clear that there is no 
significant distortion of the microwave pulse phase with 
the correction, as shown in the right panel of Supplemen¬ 
tary Figurej^. In Supplementary Figurej^, we compare 
the average gate hdelity theoretically with the microwave 
pulse with and without the phase correction. Because of 
the pulse distortions, the BBlinC 7 r /2 gate becomes less 
error-resilient, without any region of hdelity higher than 
0.9999 (see the left panel of Supplementary Figure ). 
The result in the right panel of Supplementary Figure]^ 
shows a region with hdelity higher than 0.9999 with the 
phase correction. 


Microwave amplitude correction. In the two- 
qubit experiment, microwave and radio-frequency pulses 
are combined with a diplexer (Marki DPX-I). We hnd 
that the leakage and rehection of the diplexer ports cause 
extra distortions of the microwave pulses. Supplemen¬ 
tary Figure [TO^ shows that there are distortion of the 
microwave amplitude. After inserting a 10 dB attenuator 
between the microwave components to suppress the leak¬ 
age and rehection, the waveform of the pulse sequence is 
improved to be close to the ideal case, as shown in Sup¬ 
plementary Figure [T 0 } 3 . 

The distortions of amplitude shown in Supplementary 
Figure [T 0|3 can be further corrected by pulse hxing tech¬ 
nique. Similar to that shown in Supplementary Figure 
|9ji, the distortions are recorded by an oscilloscope, and 
then fed back to the arbitrary waveform generator so that 
the distortions are minimized. Supplementary Figure [TT] 
shows the comparison of the pulse waveforms without 
and with amplitude correction. 


G. Supplementary Note 7 

Normalization of the experimental data. In the 

single-qubit experiment, the normalization is carried out 
by performing a nutation experiment |12j . The normal¬ 
ized data corresponds to the population of | 0 ) for the 
final state. 

In the two-qubit experiment, the population of \ms = 
0 ,TO/ = 1 ) (Fioi)) for the hnal state is obtained by nor¬ 
malization. According to the Ref. [T31 each occupied en¬ 
ergy level contributes to the measured photoluminescence 
intensity (/pu) with a different PL rate and these differ¬ 
ent PL rates are measured and used to determine the 
population of the levels with several sequences. Herein 
we introduce an alternative method for normalization. 
The pulse sequences for the normalization are shown in 
Supplementary Figure [T^ The measured /pp is 

Iph = ^l3\k)Pp,\k), (27) 

|fe) 

where \k) denotes the nine energy levels of the NV center 
{\ms,mj) with ms = 0 , ±1 and mj = 0 ,± 1 ), Pp^\k) is 
the population of \k) for the state p, pk is the photolu¬ 
minescence intensity if the state is \k). 

In Supplementary Figure pi denotes the initial¬ 
ized state after initializing laser pulse, pf denotes the 
final state after applying control sequence to pi. The 
RF (RFl) TT pulse exchanges the population of \ms = 
0 , TO/ = 1) and \ms = 0 , to/ = 0 ) (|tos = 0 , to/ = 0 ) and 
\ms = 0,TO/ = —1)). The measured /pp after the four 
sequences (S'/, S' 2 , i?i, and i? 2 , respectively) satisfy 

Si — S 2 = (/3|o,i) ~ /3|o.-i))(7^|o.i) ~ 7’|o.-i))) (28) 

i?i — i ?2 = (/?|o.i) ~ /3|o,-i))(7’pi.|o,i) ~ Ppi,\o,-i))j (29) 
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where P\k) (-Fpi,|fe>i) is the population of |fc) for pf (pi). 
After the initializing laser pulse, the electron spin is po¬ 
larized with coefficient a, and the nuclear spin is almost 
completely polarized. Thus we have T’pjjo,!) = ol and 
Pp,jo,-i) = 0. The state \ras = 0,m/ = —1) remains 
idle during the control sequence, thus T]o,-i) = 0- The 
population of jms = 0, mj = 1) for the final state can be 
derived as 


T’loi) = ot 


S 1 -S 2 

R 1 -R 2 ' 


(30) 


H. Supplementary Note 8 


Measuring the polarization of the NV electron 

spin. The measurement of the polarization described 
here is similar to that described before m- Supple¬ 
mentary Figure shows the results and pulse sequences 
used in the measurement. We first recorded the nuclear 
Rabi oscillation by driving the \ms = 0,mi = 1) and 
I ms = 0 ,m 7 = 0) transition. The nuclear spin is almost 
completely polarized. The amplitude of this nuclear Rabi 
oscillation is proportional to the polarization a of the 
electron spin, with 

~ /3|o,o>)«- (31) 

Secondly, another nuclear Rabi oscillation is recorded af¬ 
ter a MW2 TT pulse. The MW2 tt pulse exchanges the 
population of |ms = 0,m/ = 1) and |ms = —l,m/ = 1). 
The amplitude of this nuclear Rabi oscillation is propor¬ 
tional to the population of |ms = —l,m 7 = 1) for the 
initialized state. 


^2 — (/3|o 4) — /3|o,o>)(l ~ Q^)/2, 

Then the polarization a can be obtained with 

1 

^ ~ 1 I o ’ 


(32) 


(33) 


With the results shown in Supplementary Figure [TSl we 
estimated the polarization of the NV electron spin to be 
a = 0.83(2). 


I. Supplementary Note 9 

Measurement of single-qubit average gate fi¬ 
delity. We first describe the method for measuring the 
average gate fidelity of single-qubit gates. The average 
gate fidelity of single-qubit gates are measured with ran¬ 
domized benchmarking (RB) method [H]. Unlike that 
with quantum process tomography, the measured fidelity 
with RB method is not limited by errors in state prepa¬ 
ration and measurement. The qubit is initialized to |0), 
then a predetermined sequence of randomized computa¬ 
tional gates is applied. Each computational gate consists 


of a Pauli gate followed by a (non-Pauli) Clifford gate. 
Pauli gates are randomly chosen to rotate the qubit about 
the ±a;, ±j/, or ±z axes for an angle tt on the Bloch 
sphere, or to be a ±J identity gate; Clifford gates are 
randomly chosen to rotate about the ice or ±j/ axes for 
an angle 7r/2. The gate sequence is followed by a final 
Clifford gate chosen to ensure that the final qubit state is 
|0) if all the gates are ideal. The fidelity of the final state 
Pf, F = (0|pf|0), is measured. The measured final state 
fidelity is averaged over different randor n se quences. The 
averaged fidelity, E, is fitted with Eqn. 
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F = 1/2-k 1/2(1-dif)(l-2eg)'. 


(34) 


where I is the number of computational gates. Eg is the 
average error per gate, and dif describes errors in state 
preparation and measurement. The average gate fidelity 
is 


F = \ — r 


(35) 


In the experiment, ±a;, ±j/ rotations are realized by 
proper microwave settings, and ±z rotations are imple¬ 
mented by a rotation of the logical frame of the qubit for 
the subsequent pulses [ISHII]. 

For the naive pulse, each Clifford gate is performed by 
a rectangular tt/2 pulse and each Pauli gate by a rectan¬ 
gular TT pulse; For the five-piece SUPCODE pulse, each 
Clifford gate is performed by a five-piece SUPCODE 2.57r 
puls e (equ ivalent to 7r/2 in the single-qubit case, see Sec¬ 
tion IV C) and each Pauli gate by a pair of five-piece 


SUPCODE 2.57r pulses; For the BBl (BBlinC) pulse, 
each Clifford gate is performed by a BBl (BBlinC) tt/2 
pulse and e ach P auli gate by a BBl (BBlinC) tt pulse 
(see Section IVC I. 

The RB results for naive, five-piece SUPCODE, BBl 
and BBlinC pulses are shown in Fig. 2(b) in the 
main text and summarized in Supplementary Table [l] 
The measured average gate fidelities are 0.99968(6), 
0.99916(8), 0.999945(6) and 0.999952(6), respectively. 


Measurement of two-qubit CNOT average gate 
fidelity.In the following we describe the method for mea¬ 
suring the average gate fidelity of two-qubit CNOT gate. 
The average gate fidelity of CNOT can be measured with 
Eqn. |17| This requires the full knowledge of the quan¬ 
tum operation which is usually very difficult to be 
obtained. Quantum process tomography has been devel¬ 
oped to characterize the quantum gates. However, this 
procedure requires a number of measurements that scale 
exponentially with the number of qubits, and the mea¬ 
sured process matrix is sensitive to errors in state prepa¬ 
ration and measurement. Randomized bechmarking and 
related techniques are developed to obtain the average 
gate fidelities. However, in the hybrid system composed 
of electron and nuclear spins, single-qubit gates on the 
nuclear spin cost longer time than the electron coherence 
time. The error of gates on the nuclear spin will dominate 
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the fidelity decay in randomized benchmarking, and the 
gate fidelity of CNOT can not be precisely determined 
this way. 


Herein, we present a method to estimate the average 
fidelity of CNOT gate. We determine the fidelity by re¬ 
peated application of the CNOT gates on the system. A 
wealth of information can be obtained by studying the 
state dynamics under repeated application of quantum 
gates [11]. In Ref. [H], CNOT gates were repeatedly 
applied on the input state generated by X_t ^/2 ® I, and 
the fidelity of final states were measured. The fidelity 
Fs decays as the number of the CNOT gate, N, is in¬ 
creased. The maximum value of N was 12 in that work. 
By assuming that the decay obeys an exponential model, 
the gate fidelity Fg can be extracted. 


The pulse sequence used in our experiment is shown in 
the inset of Fig. 4e in the main text. The initial state 
of the two-qubit system is prepared by applying a RF 
tt/2 pulse after the initial laser pulse. Then N, which 
is even, times of repeated CNOT gates are applied. Fi¬ 
nally, the population of state |01) (F|oi)) after another 
RF tt/2 is measured. Up to 192 CNOT gates are ap¬ 
plied, the dynamics of however, does not obey a 

simple exponential decay. As shown in Supplementary 
Figure |T4| the measured P|oi) oscillates while decaying 


with A. ih our experiment, the nuclear spin qubit is ex¬ 
tremely ‘clean’ due to being insensitive to the external 
noises. The CNOT gate designed by quantum optimal 
control method consists of microwave pulses only. Thus 
the decay is due to the static fluctuation of Jq and (5i, 
while the oscillation is mainly due to the deviation of the 
experimental operation from the ideal CNOT gate. 


We simulated the dynamics of P|oi) based on the 
Hamiltonian Protj the pulse sequence (shown in Sup¬ 
plementary Figure and the distributions of Sq and 
(5i. The simulated dynamics of P|oi> is presented as the 
blue dashed line in Supplementary Figure |14[ The de¬ 
viation between the experimental result and the simu¬ 
lated result shown in Supplementary Figure[^ is mainly 
due to the difference between the and the prac¬ 

tical one i?rot,prac in the experiment. The microwave 
frequency does not equal the resonance frequency ex¬ 
actly, so the off-resonance term (5H is of a prior unknown 
nonzero value. The value of hyperfine coupling strength 
A = —2.16 MHz, which is used for pulse sequence design¬ 
ing, can also deviate slightly from the practical one Aexp- 
We denote this difference as SA = Aexp — A. The practi¬ 
cal Hamiltonian i?rot,prac can be extracted by fitting the 
experimental data. The fitting procedure is accomplished 
with Matlab. The best-fit values of the parameters are 
SA = 0.008(1) MHz and = 0.068(4) MHz, with the 
errors being the uncertainty within 95% confidence. The 
extracted values of SA and SO, are much smaller than the 
value of CW spectrum’s linewidth. The fitting result, 
which agrees with the experimental data well, is shown 
as the red solid line in Supplementary Figure [Tij With 
the values of SA and SO, , we can determine the fidelity of 
CNOT gate according to Eqn. [T^ The derived average 


gate fidelity is 0.9920(1), where the error is due to the 
uncertainty of SA and S^l. 


J. Supplementary Note 10 

Robust and precise optimal control method on 
NV-NV system. We have demonstrated a high fidelity 
CNOT gate at fault-tolerant threshold, taking the NV 
electron spin and ^“^N nuclear spin as qubits. The CNOT 
gate is designed with modified optimal control method. 
This method can also be used to design robust and pre¬ 
cise quantum gates on NV-NV coupled system, a key 
ingredient for scalable quantum computation using dia¬ 
mond. 

The static Hamiltonian of two coupled NV centers can 
be described as 


Hq = RnV,! + RnV, 2 + Rint, 

(36) 

^^NV,1 = 2,TtDSIi — 7eRo.l • Si, 

(37) 

H^Y ,2 = 2 ttDS12 — 'leBo,2 ■ S 2 , 

(38) 

Rint = 2 TrSi ■C-S2, 

(39) 


where Si and S 2 are the spin operators of individual NV 
centers, NV 1 and NV 2, respectively. The zero filed 
splitting is D =2870 MHz. The coupling tensor between 
NV 1 and NV 2 is denoted as C. The static magnetic 
field applied on NV 1 (NV 2) is Ron (- 80,2 )■ 

The system can be controlled by oscillating magnetic 
fields. The corresponding control Hamiltonian is 

He, NV-Nv{t) = 

- 7e ^ COs[2TTfmt + (j)m{t)]Bi^rn{t) ■ {Si + S 2 ), 
m 

where fm are the carrier frequencies of the control fields, 
Bi^rn contain the amplitudes = \Bi,m\ and the po¬ 
larization Um = Ri,m/Ri,m- The amplitudes Ri,rn and 
the phases c/m can be changed in time to steer the system. 

We turn to a rotating frame, in which the Hamiltonian 
is 

H'{t) = 

(41) 

The evolution operator with a time duration T is 

t7e(r) = (42) 

where T is the time-ordering operator. 

Similar to that described in Section |IVD[ the pulse 
sequence for a target two-qubit unitary gate U can be 
designed by maximizing the performance function 

( 43 ) 
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with d =4 and 


M' = U'frUe{T)r, (44) 

where V is projection operator on the two-qubit sub¬ 
space. 

Considering the quasi-static noises from environment 
and the control fields, the control Hamiltonian in Eqn. 
[40|is replaced by 

He, err, NV-Nv(i) = 27ri5o,l<S'zl + ‘2tt6o^2Sz2 

-7e^cos[27r (f)](l -|- ^l,m, rel) 

m 

X Bi^rn{t) • (Si + S2), 


of qubit 1 (qubit 2). Microwave pulses with two frequen¬ 
cies, which are resonant frequencies for the two qubits, 
are applied to control the system. The amplitude and 
phase of each microwave pulses can be modulated to re¬ 
alize the CNOT gate. The CNOT gate is designed so 
that the state of qubit 2 is flipped iff qubit 1 is in state 
|1). To make the CNOT gate robust against the noises, 
we consider the quasi-static distributions of the noises in 
the optimization of the pulse sequence. The performance 
function is defined as the average gate fidelity 


F' 

a 




/ 


dSl,2, rel 


With He, NV-Nv(i) replaced by He, err, NV-Nv(f), it is 
straight forward to calculate the evolution operator 
Ue, err(iC) and gate fidelity F^^^^ with quasi-static 
noises described by So,i, ^0,2) and rei- Then the 
performance function is defined by integrating g,.,, 
over distributions of ( 5 o,i, do,2, and rei- By maximiz¬ 
ing the performance function, pulse sequence for target 
F can be designed to be robust against the noises. 

We take the optimization of pulse sequence for a ro¬ 
bust and precise CNOT gate as an example. The static 
magnetic field applied on each NV center is aligned along 
the NV symmetry axis, and the NV centers can be indi¬ 
vidually addressable by application of gradient magnetic 
field. The coupling strength is taken to be 100 kHz, cor¬ 
responding to a distance of about 8 nm between two NV 
centers [20] ■ According to Ref. m, a magnetic-field 
gradient of 12 G nm“^ is available, corresponding to a 
difference of more than 200 MHz between the NV cen¬ 
ters’ resonant frequencies. The spin states \ms = 0) and 
\ms = —1) of NV 1 (NV 2) are encoded as |0) and |1) 


X Co,1(^0,1)^1.1. rel(^l,l, rel)Cb,2(<5o,2)C’l,2, rel(<5l,2, rel) 

X F' 

seq, err ’ 

(46) 

where Po.i(<^o,i) and Pyi, rei(^i,i, rei) {Po,2{So,2) and 
Pi, 2 , rei(<5i,2, rel)) desci'ibe the distributions of quasi-static 
dephasing noise and control field fluctuation on NV 1 
(NV 2). The distribution Co. 1(^0,1) (Cb, 2 (^o. 2 )) is mainly 
due to the ^^C nuclear spin bath surrounding the NV 1 
(NV 2). Since is naturally abundant, the distribu¬ 
tion Po{6o) in the main text is a typical estimation of 
Cb,i((5o,i) iPo,2{5o,2))- The distributions rei(i5i.i. rei) 
and Pi^2, rei(i^i,2, rel), mainly depending on the microwave 
generator, are consistent with Pi(i5i) which can be ob¬ 
tained from experiment data. Considering the distribu¬ 
tions, a pulse sequence for the CNOT gate can be op¬ 
timized to achieve an average gate fidelity F^ = 0.9926 
by our method. The pulse sequence is shown in Supple¬ 
mentary Figure [15| Thus our method can be applied to 
realize robust and high fidelity two-qubit gate on spatially 
separated NV centers. 
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